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ABSTRACT 


The  Green’s  function  approach  for  the  subsonic  compressible  potential 
flow  over  a hel< copter  rotor  in  hover  or  steady  state  vertical  climb  was 
examined.  Several  methods  of  approach  were  taken  to  derive  the  Green's 
function.  An  attempt  to  derive  the  Green’s  function  directly  from  the 
differential  equation  written  in  the  non-inertial  coordinate  system 
(rotating  Cartesian  coordinates  attached  to  the  blade)  was  not  successful. 
However,  since  the  Greenfe  function  was  known  in  the  inertial  system,  a 
transformation  of  the  integral  equation  for  the  velocity  potential  in  this 
system  was  transformed  successfully  to  the  blade  fixed  coordinate  system. 
The  resulting  integral  equation  involves  surface  derivatives  of  (p  which 
precludes  the  use  of  previous  numerical  :>chemes  for  its  solution. 


A second  part  of  the  report  includes  the  results  of  the  second  iteration 
for  the  force-free  wake  of  a rotor  in  hover. 
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1.  INTRODUCTION 


"■*  daman  1 for  heavily  loaded,  high  performance  helicopters  has 
/'•w# r. «4i  # great  deal  of  research  work,  on  the  aerodynamics  of  the 
• •*  Lifting  line  theory  has  been  employed  to  Indicate  the 

•fleets  ef  the  distorted  wake  vortices.  Various  ways  of  prescribing 
-*  t>;«ttt*  the  distorted  wake  were  developed,  however,  due  to  the 
i •weemtu  .a  MMt  of  calculation,  very  few  studies  have  been  made 
.4,4  i • ' tag  surface  theory  applied  to  helicopter  rotor  blades. 

M] at  drawback  of  a lifting  line  model  is  that  it  cannot 
pnifeac  lit  variation  of  the  Induced  flow  In  chordwlse  direction  which 
■«  <**>  • tan  If leant  for  a heavily  loaded  helicopter,  especially  in 

.if ting  surface  theory  can  give  the  details  of  the  flow  in  both 
.......  .=  ».- ( -.,r  Ilrfrtlors.  This  is  extremely  important  in  new 

Mill  dtoalga  hoc sues  the  method  la  capable  of  showing  the  effects  of  any 
■hawse  la  Slade  geoaetry. 

1MN  arm  three  methods  of  approach  that  can  be  used  when  applying 
:«'*(**  surface  theory.  One  method  of  approach  predicts  the  local  surface 
If  IMSdag  a loading  function  which  Is  expressed  in  a series 
af  aaSMSMl  madam  with  walmmm  coefficients.  These  unknown  coefficients 
Spa  Urnm  attained  by  aatlefylng  the  normal  velocity  conditions,  either 
IMntCJy  rn  Imdlrectly  at  a eat  of  points  on  the  surface.  The  second 
ea»a*«t  »f  apyreeeh  predicts  the  local  surface  loadings  by  first  assuming 
sw  Isfleamea  tenet Ion,  such  as  sink  and  source,  vortex,  or  doublets,  of 
igiln  distributed  on  the  surface.  By  satisfying  the  normal 
JMaett?  eaadltiea,  so  integral  equation  is  formed  and  then  solved  by 
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approximating  the  surface  with  a finite  number  of  elements.  The  solution 
of  the  influence  distribution  provides  the  velocities,  pressures  and 
hence  the  loadings.  The  third  method  predicts  the  local  surface  loading 
by  directly  solving  the  governing  elation  in  terms  of  the  velocity 
potential  through  the  application  of  Green's  theorem.  An  integral  equation 
is  formed  and  solved  in  the  same  way  as  in  the  second  method.  In  fact, 
this  method  results  in  the  same  integral  equation  as  the  second  method 
when  the  flow  is  incompressible.  However,  the  use  of  Green's  theorem 
generalizes  the  second  method  and  is  capable  of  solving  compressible 
flow  problems. 

It  should  be  noticed  that  the  first  method  gives  the  pressure 
coefficient  without  going  through  the  calculation  of  velocities  and 
pressures,  as  is  required  in  the  latter  two  methods.  When  the  actual 
wetted  surface  of  the  blade  is  used  to  satisfy  the  normal  boundary 
conditions  and  to  calculate  the  velocities  and  pressures,  these  data 
can  be  used  in  the  analysis  of  viscous  flow  problems. 

An  attempt  has  been  made  to  apply  lifting  surface  theory  using  the 

second  approach  to  an  arbitrary  helicopter  rotor  blade  system  which  is 

in  axial  motion  including  hover  as  a special  case.  A test  case  was  run 

with  the  elements  assumed  on  the  actual  wetted  surface  of  the  blade. 

Two  reports  have  been  published  on  this  work.  The  first  report^  assumed 

a classical  helical  rigid  wake  with  a doublet  distribution  on  the  wetted 

2 

surface  of  the  blade  and  the  wake.  The  second  report  attempted  to 
compute  the  force  free  wake  by  starting  with  the  rigid  wake  results  and 
iterating  on  the  geometry  of  the  wake.  A force  free  wake  requires  that 
each  point  on  the  wake  is  free  to  convect  at  its  local  velocity.  However, 
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because  of  the  complexity  of  this  problem,  only  the  results  of  the  first 
iteration  were  included. 

Both  of  the  above  referenced  reports  are  concerned  with  incompressible 
flows.  In  order  to  investigate  the  effect  of  compressibility,  in  this 
report  the  compressible  case  of  helicopter  rotor  lifting  surface  theory 
will  be  considered.  Since  the  influence  functions  for  an  incompressible 
case  will  not  satisfy  the  governing  equation  of  a compressible  case,  the 
third  approach  was  used.  In  other  words,  Green's  theorem  was  applied  in 
this  case  in  order  to  get  an  appropriate  Influence  function.  The  analysis 
will  be  compared  with  a general  compressible  theory.  In  addition,  the 
iteration  scheme  used  in  Ref.  2 was  further  pursued  and  modified,  and 
the  results  of  the  second  iteration  are  included  in  this  report. 
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2.  LITERATURE  REVIEW 


Extensive  reviews,  as  well  as  listing  of  references  on  lifting 
surface  theory  and  helicopter  wake  analysis  were  presented  In  Refs.  1 
and  2.  In  summary,  very  few  reports  have  been  published  on  the 
application  of  lifting  surface  theory  to  helicopter  rotor  blades.  Most 
of  the  previous  studies  circumvent  the  complicated  calculations  by  using 
the  camber  plane  Instead  of  the  wetted  surface  of  the  blade.  The 
description  of  helicopter  wake  geometry,  either  the  prescribed  rigid 
wake  or  force  free  distorted  wake,  has  had  some  success  under  the 
limitations  of  the  lifting  line  model.  In  most  cases  the  wake  sheet 
concept  were  not  maintained  wherein  the  tip  vortex  is  separated  from 
the  Inboard  vortices.  In  the  compressible  helicopter  rotor  lifting 
surface  theory  to  be  considered  in  this  report,  the  true  wetted  surface 
on  the  blade  will  be  used  and  the  sheet  concept  will  be  maintained.  The 
Green' 8 theorem  approach  will  be  used  to  get  the  influence  function. 

The  literature  concerned  with  the  compressible  flow  over  a helicopter 
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rotor  blade  system  is  sparse.  In  1969,  Sopher  presented  an  analysis  on 

a hovering,  nonlifting  thickness  problem.  The  solution  of  the  velocity 

potential  in  the  general  aeoustic  equation  was  taken  from  Garrick's  non- 
4 

steady  wing  theory.  With  some  approximations  made  for  the  steady  case, 
the  velocities  and  pressures  were  calculated  and  then  compared  with  the 
results  of  blade  element  theory.  It  illustrated  that  blade  element 
theory  was  inaccurate  near  the  blade  tip  and  became  worde  as  the  blade 
speed  was  increased.  This  was  the  first  attempt  to  check  the  dependence 
of  three-dimensional  flow  on  compressibility. 
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Caradonna  and  Isom'’  also  solved  a hovering,  nonlifting  problem  by 
using  the  mixed-difference  relaxation  method.  Both  sub3onlc  and  transonic 
flows  were  considered.  Scale  factors  in  terms  of  aspect  ratio  were 
introduced  to  simplify  the  equations.  These  equations  were  then  solved 
by  a finite  difference  scheme.  The  results  show  the  importance  of  tip 
Mach  number  and  aspect  ratio  on  the  growth  and  extent  of  shock  waves  in 
the  tip  region,  add  indicate  a significant  reduction  in  shock  strength 
with  decreasing  aspect  ratio. 

A second  paper**  was  presented  by  Caradonna  and  Isom  in  1975  on  the 
transonic  case  in  formard  motion.  Scale  factors  were  employed  and 
Successive  Over  Relaxation  (S.O.R.)  schemes  were  used.  This  was  an 
unsteady  problem.  The  results  show  that  the  flow  in  the  tip  region  is 
mostly  unsteady  in  the  decelerating  flow  region.  The  influence  of  aspect 
ratio,  advance  ratio  and  Mach  number  on  this  process  was  discussed. 

All  of  the  papers  cited  above  (Refs.  4,  5,  and  6)  deal  with  nonlifting 
problems  and  these  are  the  only  major  works  that  have  attempted  to  solve 
the  compressible  flow  over  helicopter  rotor  blades.  The  approaches  are 
different  from  what  is  ususally  called  lifting  surface  theory  and  wakes 
are  not  included  because  there  is  no  lift. 

A general  formulation  of  subsonic  lifting  surface  theory  using  the 
Green's  theorem  approach  was  given  by  Morino  and  Kuo.^  The  theory  is  so 
general  that  it  can  include  any  type  of  motion,  even  rotational.  However, 
it  is  very  difficult  to  visualize  the  rotational  problem  unless  the 
coordinate  system  is  attached  to  the  rotating  body.  No  attempt  was  made 
to  apply  this  theory  to  helicopters.  In  this  report,  this  general  equation 
will  be  transformed  into  rotational  coordinate  systems  in  order  to  compare 
with  the  present  analysis. 
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3.  FORMULATION  OF  THE  PROBLEM 


3.1  Problem  Definition 

Steddy  potential  flow  Is  considered  over  an  arbitrary  helioopter  rotor 
with  any  number  of  blades  In  axial  motion  Including  hover.  The  rotational 
speed  of  the  rotor  is  assumed  to  be  large  so  that  compressibility  effects 
must  be  taken  into  account.  The  problem  is  to  determine  the  velocities 
and  pressures  on  the  wetted  surfaces  of  the  blade  and  to  study  the  effects 
of  compressibility.  The  influence  of  the  wake  will  also  be  Included. 


3 . ? Governing  Equation 

If  a Cartesian  coordinate  system  is  attached  to  one  of  the  rotating 
blades,  (see  Fig.  1),  the  governing  equation  for  linearized  subsonic  or 

l 

supersonic  flow  was  shown  by  Caradonna  and  Isom'’  to  be  as  follows: 


n * »i 


(3-1) 


where  CL*o  is  the  sonic  speed  of  the  medium  at  rest, 
lO  is  the  rotational  speed  of  the  blade,  and 
f is  the  perturbation  velocity  potential,  defined  in 
the  following  equation  . 

The  total  velocity  at  any  point  ia  the  flow  field  is 


+ Vp  f 


(3-2) 
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■** 


i! 


-♦  4 A « 

where  Y=  XX+tyy  + 2%,  Yff  means  taking  gradient  with  respect  to  the  XfZ 

coordinate  system,  and  1 ^ is  the  climbing  speed.  For  the  hover  case 

- 0.  V is  the  velocity  as  viewed  by  an  observer  in  the  rotating 
coordinate  system. 

3 . 3 Boundary  Conditions 

The  boundary  conditions  are  the  same  as  in  the  incompressible  case. 
That  is, 

(1)  At  very  far  regions  where  there  is  no  induced  velocity 

Vr?**0-  o-3) 

(2)  On  the  blade  surface,  the  normal  velocity  is  0 

7*  = v-  n.  - 0 . (3-4) 


3.4  Pressure  Coefficient 

The  pressure  coefficient  is  usually  defined  by 

P-K 
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(3-5) 


where  from  Eq.  (3-2)  is  equal  to  \ - (uJ £ X T ) + )/x  $ | , which  is 

dependent  on  the  radius  Y ■ In  order  to  use  a constant  as  the  reference, 
Cr  is  redefined  as 

„ P-K 


T V*) 


(3-6) 


where  R is  the  rotor  tip  radius.  For  the  hover  case 

t-K 


C.=- 


If  the  parameter  0 is  defined  as 

u > R 


(3-7) 


?- 


(3-8) 


By  using  perfect  gas  relations 


To  express  j>  in  terms  of  the  velocity,  one  has  to  consult  the  energy 
equation.  For  compressible  or  incompressible  flow 


where  Cp  is  the  specific  heat  and  subscript  0 designates  the  reservoir 
condition.  With  the  expressions  Ct  /(¥’’“ /)  and  = i* RT 


This  can  be  rewritten  as 
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Dividing  Eq.  (3-15)  by  Eq.  (3-14)  given 


t _r  * + irt)  l'-' 
& L z + (r~i)( _v_ J 


(3-16) 


To  eliminate  , one  may  write  Eq.  (3-12)  in  another  form 


v/  a 

— 1-  

-4  /'-I 


1 V* 
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j ^Ar 


r~\ 


(3-17) 


This  gives 


Y-\ 


r-\ 


(3-18) 


or 


/ /»  - 1 - V — 

• 1 


2 v ■>  i1 

Substituting  this  into  Eq.  (3-16)  one  gets 


(3-19) 


(3-20) 


therefore  the  pressure  coefficient  from  Eq.  (3-10)  is 


L-i 

fl+  2 

pk  f 

f * 

11  iC'J  '/ 

(3-21) 


This  is  the  exact  expression  for  . Since  the  local  velocity 
can  be  easily  obtained  from  Eq.  (3-2),  there  is  no  need  to  expand  this 
out  to  get  an  approximation.  It  should  be  noted  that  for  a rotor  blade 
usually  has  a significant  magnitude  in  the  spanwise  direction 


as 
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well  as  In  the  chordwise  direction,  and  is  a variable  depending  on  the 
radius.  It  is  dangerous  to  use  the  familiar  expression  for 

regular  wings  especially  near  the  inboard  tip  because  the  contribution 
from  the  spanwise  direction  will  be  overlooked. 


3.5  Nondimens ionalization 

It  is  always  convenient  to  nondimensionalize  the  equations  and 
introduce  some  nondimensional  parameters,  as  given  by  Eq.  (3-8) 

is  a very  important  parameter  for  compressible  rotary  flow  and  is 
equivalent  to  the  Mach  number  in  rectilinear  flow.  It  will  be  used 
throughout  these  equations. 

In  the  following,  all  the  lengths  will  be  divided  by  tip  radius  ft 
end  all  the  velocities  will  be  divided  by  tip  speed  u)R,  Let  superscript 
' designate  nondimensional  quantities,  and  define 


n 


( 


X 

R 

JL. 

R 

2 

R 

UJ/f 

n 

R 

jL 

Jk. 

UfR 


etc. 


J 


(3-22) 


The  governing  equation,  Eq.  (3-1)  can  be  rewritten  as 
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The  boundary  conditions  become 
(1)  in  free  stream 


V?  ?'=  o 

(2)  on  blade  surface 


(3-23) 

(3-24) 

(3-25) 

(3-25a) 


Other  quantities  and  equations  can  be  written  in  a similar  way. 
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4.  APPLICATION  OF  GREEN’S  THEOREM 


For  the  sake  of  convenience,  the  superscript  "/"  in  the  ncndimensional 
equations  will  be  dropped  in  this  section,  without  being  confused  with 
the  dimensional  ones.  Eq.  (3-23)  may  be  rewritten  as 


-+  %- 


= L f 

The  boundary  conditions  now  have  the  form 
(1)  in  the  free  stream 


d X 


(4-1) 


VrT  <P  = 0 (4-2) 

r > 

(2)  on  blade  surface 


To  solve  Eq.  (4-1),  one  may  multiply  the  equation  by  a function  ^ , 
which  is  called  Green's  function  and  will  be  determined  later,  and 
integrate  it  throughout  the  volume  in  which  the  flow  field  is  being 
considered,  then  transform  it  into  an  integral  equation.  The  volume 
of  integration  contains  the  region  from  the  body  surface  to  infinity, 
therefore  this  is  an  exterior  problem.  Integration  by  parts  was  applied 
to  each  term  and  the  volume  integral  was  transformed  into  surface  integral. 
This  is  actually  the  approach  of  Green's  theorem.  The  details  of  this 
derivation  are  given  In  Appendix  A. 

Fig.  2 shows  the  volume  of  integration  y and  the  surface  boundary 
5 . The  result  of  Green's  theorem  is 
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JjjfrLf  dr 

+ { *[  i - ( p <1  *"b-u*r]ift-*t(*)  **&VM » 
+ { 4-lf  ~tf  ?}  ;J  • J 


*=  0. 


(4-4) 


Let  the  operator  in  the  curled  bracket  of  the  volume  integral  be 

L*  - [ i- - <»?]&  f [i- lf  + 7F 


Incidentally,  by  rearranging,  one  finds 


(4-6) 


which  turns  out  to  be  the  same  as  |_  , i.e.,J_  =■]_  • This  shows  that 
L is  self-adjoint. 

If  one  can  find  the  function  such  that 

L*<k-L$-  S(z-x,,  <4-7> 

where  £( X~~*h  fli  > i~ 2t)  is  the  Tirac  delta  function,  then  by 


definition  of  Dirac  delta  function,  the  volume  integral  becomes 


(4-9) 

This  includes  the  surface  integral  only.  It  could  be  solved  by  approxi 


mating  the  surface  with  a finite  number  of  elements  if  Green's  function 
Q-  can  be  obtained.  The  remaining  question  is  how  to  find  £ . 

Gr  has  to  be  solved  from  Eq.  (4-7).  The  symmetry  of  the  operator 


suggests  that  it  is  better  to  do  this  in  cyclindrical  coordinates 


(4-10) 

Consider  first  and  separate  the  variables.  The  details  are 

shown  in  Appendix  B.  It  turns  out  that  /_  is  of  the  standard  Bessel  type 


and  Qf  may  be  expressed  in  terms  of  Bessel  functions.  Unfortunately, 
it  still  has  to  be  written  in  integral  form  and  no  compact  expressions 


can  be  achieved 


Several  attempts  were  made  to  expand  the  integral  form  in  series 


such  as 


however  all  these  series  are  divergent 


The  technique  of  perturbation  was  also  applied  in  light  of  the  fact  that 


r 

L 

L 


15 


t In*  toe  (he  Laplace  equation  is  well  known,  but  singularity 
iietuw  «»,i h.*a( •( *4 . TLe  details  of  all  of  the  above  investigation 

toto  000  a*  tHMMti  (a  (his  report  because  of  its  complexity. 

fh*  *«M4  lee  Opean's  function  has  been  one  of  the  main  tasks  of 
toto  PMfMI*  iOMpeef • e werbeble  tore  could  not  be  obtained  and 
t>t  |*  »)  U *M  eeleeble  te  Its  present  fore. 

» an**>»U  be  neat  loeed  finally  with  regard  to  notation, 
to  DM*  MCCtoo  lie  deear  portables  of  integration  do  not  have  the 
Mtototto  “I',  eftoteae  tie  coordinates  of  the  field  point  do.  This 
to  ♦ MMta  ed  the  derteotlee  of  Orson* a theoren,  and  is  opposite 


5.  COMPARISON  WITH  A GENERAL  COMPRESSIBLE  THEORY 


Unsteady  compressible  flow  was  treated  by  Marino  and  Kuo  by  using 


Green’s  function.  This  is  a very  general  theory  that  will  accomodate 


any  complex  configuration  in  any  motion.  However,  in  order  to  visualize 


the  rotational  motion  of  the  blades,  one  has  to  transform  the  equation 


This  will  be  done  in  this  section 


The  perturbation  velocity  potential  in  Marino  and  Kuo's  paper 


(designated  as  MK  in  the  following)  is  different  from  that  in  this 


which  is  the  free  stream  velocity  of  the 


body.  If  the  total  velocity  potential  is 


the  perturbation 


potential 


is  defined  by 


while  this  report,  from  Eq.  (3-2),  uses 


Clarifying  this  difference,  one  can  rewrite  the  governing  equation 
in  MK.  t with  respect  to  an  Inertial  coordinate  system  Xy  € £ 
attached  to  the  body  and  moving  with  velocity  (J«,  In  x-direction,  as 


where  all  the  notations  have  been,  and  will  be,  written  in  this  report’ 
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convention.  The  surface  of  integration  is  represented  by 

SC*,  a,  *)=*0  (5-5) 

The  boundary  condition  can  be  written  as  P5 /P£  =*  0 , or 


i_ 

>»  )t^5( 


(5-6) 


All  the  overhead  bars  (— ) mean  the  quantity  in  the  inertial  coordinate 
system. 

For  the  hovering  helicopter  case,  the  inertial  coordinate  system  is 
attached  to  the  body  of  the  helicopter  with  z-axis  at  the  center  of  the 
hub.  Since  there  is  no  advancing  speed,  U^—0  » the  transformation 

to  the  rotational  X ^ systems  (refer  to 

Fig.  3)  then  has  the  following  relationship: 


from 


V 


«§T 


system 


* = X ces(u>t)+  jj  sin(<4t  ) 
y = -X  sin  (ts>£)  ■+}  cos  Cut) 

A 

2 = •? 


(5-7) 


t = >t  J 

The  time  \ in  rotational  coordinate  systems  will  be  considered  for 
the  transformation  and  ^ terms  will  be  dropped  later  for  the  steady 
case.  Using  Eq.  (5-7),  the  following  relevant  relations  can  be  obtained: 


% = X COS(u)t)  - y Sin  (U)A) 
J = * sin  (u)t)  -\-ij  coS  C u)t) 
2=2 


(5-8) 


a = r 


r w|i, 


- (tv£  *7)-  vr  -+- 


X cos(u)t)-hj  J,\  (wit) 
-"X  Sin(ivfc)  cos  (oJZ) 


X cos  (cut)  - jj  smite X) 
% Sin  (cot)  -f-  y cos  (cot) 
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The  Jacobian  of  the  transformation 


j ^ = j 


(5-14) 


Using  the  above  relations,  with  U^  — D and  =1  0 , one  can  rewrite 

dX 


Eq.  (5-4)  as 


V?r--L[-(ui  X r).  |7r]^  =0. 


(5-15) 


Rembering  (u>!X  ? ) ■ ^ ftjx  and  expanding  the  square  brackets 


in  the  above,  this  now  has  the  form 


r '-<?;]£. r.-(£rj$*£ 


(5-16) 


which  is  just  the  governing  equation  given  in  Section  3.  The  boundary 


condition  becomes 


[-(*?'?)]  -VS'  (5-1,, 

The  Breen’s  function  1h  the  inertial  system  for  the  unsteady  case  is 


Gr-~ 


£ (*,  -t  ■+  o) 


(5-18) 


where 


Cz-t'f])*  (5-i9) 


m 


In  the  rotational  system  now,  with  = 0 » p(=]  , 


4ir  I* 


but  now 


(5-22) 


Ki  — [ xi)  ~t  (<f~  'jt)  ■+  ( 4 ~~  &i)  J 


(5-23) 


When  Green's  theorem  is  applied  to  Eq.  (5-4),  the  integral  equation 
in  the  general  form  with  respect  to  the  inertial  coordinate  system  is 

inEf  (x,  J,  *,  f) 

~ fT  rv>s-Y<p--i-  — 

Jj  L 7 Vl  / nl  Jt  v . 


I 
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£=  I 


Inside  tf~ 


on  (T 


outside  (T 


(7  d 

and  [ J indicates  evaluation  atf(=f-^.  Similarly,  ^ is  the  surface 
given  by 

S0~  SlXi.JuhJ-*)  (5-27) 

where  (%fl  ljt , 2-,  ) is  the  dummy  point  of  integration  on  . 

Substituting  in  the  transformation  relations,  Eq.  (5-24)  in  the 
inertial  system  with  and  — 0 becomes 


=-Sfl  »5r  s- % 

+ If  { 7r,  $ ' vf,  (-g ) + l(  w % * r’ ) ' sJf ( ui'y r‘  Wj^j 


where  (5_28) 

(5-29) 

9 0 a* 

and  (T  or  ^ now  is  just  the  surface  given  by  $(XM,i}=0,  which  is 
independent  of  time. 

Substituting  the  boundary  condition,  Eq.  (5-17),  and  remembering 

(U/^O-Pir-  .„d  ^4-^-^:. 

''  |&J  r ZK, 

Eq.  (5-28)  can  be  rewritten  as  1 


+rEf(*.f.z) 


All  of  the  above  equations  in  this  section  are  written  in  dimensional 
form.  In  order  to  compare  with  the  results  in  the  last  section,  these 
equations  have  to  be  nondimensionalized.  Using  the  same  notation  defined 
in  Section  3,  and  noting  that  — <(x  4 J — R*  ■Ax'dy'**  ftJs'o ne  can 
rewrite  Eq.  (5-30)  as 

WE  fix',  ft) 


To  compare  this  with  Eq.  (4-9),  one  should  ignore  the  primes  ('V") 
and  remember  that  the  variables  with  subscript  are  dummies  of 
integration.  Eq.  (4-9)  can  be  rewritten  in  the  following  form  using  the 
convention  of  this  section 


r 
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f'Wf?) 


Though  this  looks  very  similar  in  form  as  Eq.  (5-31),  it  is  not  easy 
to  get  an  expression  for  Qr  by  comparison. 

The  terms  in  Eq.  (5-31)  can  be  substituted  by  the  boundary 

condition  shown  in  Eq.  (3-25a).  However,  the  term  ( — 

*{>*'/*}  ' ) is  not  known.  This  makes  Eq.  (5-31)  difficult  to  be 

solved  by  transforming  it  into  an  algebraic  equation.  This  is  also  the 
reason  why  no  numerical  solution  is  attempted  in  this  report. 
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6.  SOME  RESULTS  FROM  SECOND  ITERATION 

This  section  will  present  some  result?  for  the  second  iteration  of 
the  force-free  wake  and  a continuation  of  the  first  iteration  results  of 
Ref.  2. 

The  damped  first  iteration  wake  geometry  calculated  in  this  reference 
was  taken  to  be  the  basis  of  the  second  iteration.  In  order  to  get  a 
better  understanding  of  the  wake  distortion  near  the  trailing  edge, 
additional  steps  were  used  in  this  region.  Moreover,  it  was  found  that 
since  the  control  points  on  the  radial  lines  were  no  longer  linear  after 
distortion,  some  irregularities  would  occur  due  to  the  differentiation 
of  in  the  radial  direction.  Therefore,  before  starting  the  second 
iteration,  the  distorted  radial  lines  should  be  projected  onto  each 
radial  plane,  which  was  done  by  using  the  projection  routine  discussed 
in  Appendix  B of  Ref.  2,  to  make  them  linear  again.  The  final  input  wake 
geometry  of  the  second  iteration  looks  like  the  one  shown  in  Fig,  4.  The 
sizes  of  the  inserted  steps  are  shown  in  Fig.  6.  The  step  size  for  the 
far  wake  calculations  was  similar  to  that  used  in  Ref.  2. 

This  second  iteration  process  has  been  carried  out,  including  finding 
a new  wake,  solving  the  equations  and  finally  recalculating  the  loadings. 
In  the  calculation  of  the  new  wake  geometry,  the  contribution  to  the 
velocity  due  to  the  derivative  of  the  velocity  potential  in  the  streamwise 
direction  was  also  included.  It  should  be  understood  that  the  value  of  J> 
along  any  given  streamline  leaving  the  trailing  edge  would  remain  constant 
and  equal  to  the  wa^e  geometry  waa  correct.  This  contribution 

was  neglected  in  the  first  iteration.  All  the  techniques  developed  for 
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the  first  iteration  were  used  and  therefore  the  program  is  essentially 
the  same. 

The  resulting  geometry  of  the  second  iteration  is  shown  in  Figs.  5 
and  6.  Instead  of  rolling  up,  the  points  near  both  tips  were  washed 
downward.  The  downward  distortion  increased  the  loading  on  the  blade 
and  therefore  made  the  result  more  divergent. 

The  damping  technique  was  used  again.  The  geometry  used  to  calculate 
the  new  loading  is  the  average  of  the  one  shown  in  Fig.  5 and  the  damped 
first  iteration.  Fig.  4. 

The  A/UfE  and  loading  curve,  are  shown  on  Figs.  7 and  8 respectively. 
The  overall  thrust  coefficient  is  .005424. 

It  is  obvious  that  the  present  result  is  caused  by  the  fact  that 
the  wake  distortion  scheme  was  not  working  well  near  the  outer  edges 
of  the  wake  sheet  so  that  the  expected  tip  roll-ups  were  not  generated. 
Before  carrying  out  additional  iterations,  the  wake  distortion  calculation 
scheme  must  be  thoroughly  investigated  in  this  region. 
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7.  CONCLUSIONS 

One  of  the  objectives  of  this  work  was  to  adapt  our  lifting  surface  theory 
programs,  as  described  in  Refs.  1 and  2,  to  the  subsonic  compressible  uce. 

One  procedure  that  was  attempted  was  the  use  of  the  Green's  theorem  approach 
whereby  the  governing  integral  equation  was  transformed  into  a set  of 
algebraic  equations.  It  was  hoped  that  this  approach  would  lead  to  an 
analytical  and  exact  form  of  the  Green's  function  so  that  the  resulting 
integral  equation  could  be  solved  by  applying  the  same  numerical  scheme 
as  developed  in  Refs.  1 and  2.  The  analysis  was  not  restricted  to  any 
particular  compressible  flow  region  on  the  blade  so  that  the  effect  of 
compressibility  would  be  accounted  for  throughout  the  complete  flow  field. 

The  attempt  to  derive  the  Green's  function  directly  from  the  differential 
equation  written  in  the  non-inertial  coordinate  system  (rotating  Cartesian 
coordinates  attached  to  the  blade)  was  not  successful.  However,  since  the 
Green's  function  was  known  in  the  inertial  system,  a transformation  of  the 
Integral  equation  for  the  velocity  potential  in  this  system  was  transformed 
successfully  to  the  blade  fixed  coordinate  system.  Investigation  of  the 
resulting  integral  equation  (Eq.  0-31))  shows  that  the  appearance  of  terms 
involving  surface  derivatives  of  cj>  precludes  the  use  of  our  previous 
numerical  scheme.  Further  attempts  were  made  to  investigate  the  Green's 
function  in  a rotating  cylindrical  coordinate  system.  It  was  found  that 
the  Green's  function  for  the  appropriate  differential  equation  would  be  in 
the  form  of  an  integral  which  would  result  in  a numerical  scheme  that  was  too 
complicated  to  be  practical. 

A second  part  of  this  study  was  related  to  improving  the  force  free  wake 
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geometry  by  extending  further  the  procedures  used  in  Ref.  2 for  the  first 
iteration  wake.  Extensive  calculations  were  carried  out  in  attempting  to 
find  a convergent  procedure  which  would  result  in  the  correct  force  free 
wake  geometry.  The  results  appear  to  show  that  the  inboard  portion  of  the 
wake  was  convergent,  however,  the  outboard  portion  of  the  wake  diverged.  It 
was  felt  that  v,'is  divergent  tendency  was  caused  by  inaccurate  calculations 
of  the  induced  velocities  near  the  edges  of  the  wake  sheet.  Further  attempts 
should  be  made  at  improving  numerical  procedures  to  compute  these  induced 
velocities  in  the  region  of  the  free  edges  of  the  wake. 
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APPENDIX  A DERIVATION  OF  GREEN'S  THEOREM 


In  this  appendix,  Eq.  (4-1)  multiplied  by  function  (r 

+ £['-(  K]  ■ 

+ p‘$(2*yj&  + y-l£-)=o  (A-i) 

will  be  integrated  throughout  on  volume  , and  by  using  integration  by  parts 
the  volume  integral  will  be  transferred  into  surface  Integral. 


Consider  each  term  one  by  one.  For  the  first  term,  let 

JM'-fW’J  itdJ 


(A- 2) 


V V 

Referring  to  Fig.  2,  consider  the  integration  in  the  x-direction  first.  For 
a strip  that  does  not  pass  through  the  internal  boundary  surface  5 • the 
limit  is  from  — Oo  to  oo  , therefore  for  x-direction 

= o-  r 

l„  »x 

T>\  >\  J 7 


'~PC  ~l)6 


(A-3) 


For  a strip  that  passes  through  ^ , it  is  split  into  two  parts: 


^i(^) 


+r 

■J*  / m f)  ' a* 


X(y) 


^ -A3  ?X  J L 2X  *X 


39 


■■ci-wc-r*'*) 


.rzM)* 00 

7 „ < 


-r 
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Now  consider  the  integration  in  y-  and  z-directions.  Since 


(A-4) 


dy  d?  s*  % ’ft  ds 


the  bracketed  terms  in  Eq.  (A-4),  after  putting  back  into  the  volume  integral. 


become 


The  other  terms  in  Eq.  (A-4)  together  with  the  term  in  Eq.  (A-3)  gives, 
for  the  total  volume  y, 

fff  J-rr 


JIf  f « 


therefore 


(A— 5 ) 


By  substitution  of  Eq.  (A-2) 


IIK'-wJ-If  ^ 

r- 


+ Ml'-(PP]^dv 

V 


Similarly, 


-h-  rO  - (f*n  4f-r}  r^s 

t*r 

TT  ^ 


JfJ  6-  # MK-  *!t+£r]  * « 

r s V 


ill  rfMJJfrC-^  rf*Vs , 

in  x-  and  y-direction  it  could  be  either 
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Together  with  the  integration  in  z-direction  by  noting 


APPENDIX  B TRANSFORMATION  OF  THE  GOVERNING 


EQUATION  INTO  BESSEL  EQUATION 


Consider  the  following  equation 
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then  Eq.  (B-l)  becomes 


Again  let 


where 
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Dividing  through  by  ft&Zj 


This  is  the  standard  form  of  Bessel  equation  whose  solution  is 


